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Abstract 

The paper discusses the use of amplitude equations to describe the spatio-temporal 
dynamics of a chemical reaction-diffusion system based on an Oregonator model of the 
Belousov-Zhabotinsky reaction. Sufficiently close to a supercritical Hopf bifurcation the 
reaction-diffusion equation can be approximated by a complex Ginzburg-Landau equation 
with parameters determined by the original equation at the point of operation considered. 
We illustrate the validity of this reduction by comparing numerical spiral wave solutions to 
the Oregonator reaction-diffusion equation with the corresponding solutions to the complex 
Ginzburg-Landau equation at finite distances from the bifurcation point. We also compare 
the solutions at a bifurcation point where the systems develop spatio-temporal chaos. We 
show that the complex Ginzburg-Landau equation represents the dynamical behavior of the 
reaction-diffusion equation remarkably well sufficiently far from the bifurcation point for 
experimental applications to be feasible. 



To appear in Int. J. Bifurcation and Chaos, 7(9), 1997. 



1 Introduction 



During the last two decades, the study of spiral waves and spatio-temporal chaos in physical, 
chemical, and biological systems has received much attention. In physics, spiral patterns have 
been studied extensively in hydrodynamic syste ms ||Cross fc Hohenberg 1993 1 . Recently, spiral 



formation has also been studied in liquid crystals HCoullet et al. 1994 1. Biological examples include 



the growth of the slime mold Dictyostelium discoideum [ Newell 1983 ] and potential differences 



occuring on the surface of rabbit heart muscles [ Allesie et al. 1977 1. In chemical systems, spiral 



formation and related phenomena like target wave propagation often occur in oscillatory chemical 
reactions which take place in unstirred media (e.g. a petri dish). In this case the spirals appear 
as concentration variations over the spatial domain of the system. A well-known example of 
a chemical reaction exhibiting such phenomena, is the Belousov-Zhabotinsky reaction, whose 
properties in spatially distributed media has been discussed extensively in the chemical literature. 



For further reference, see Zhabotinsky [1991] and references therein. Recently, observations of 
spiral wave development and spatio-temporally chaotic patterns close to a supercritical Hopf 
bifurcation was reported by Ouyang & Flesselles |1996 |. 



The state space of chemical systems exhibiting wave and spatio-temporally chaotic phenomena 
almost always have a fairly high dimension. Consequently, realistic modeling of such phenomena 
by direct integration of the reaction-diffusion equation is difficult or sometimes even impossible. 
For this reason most works in this field have used simplified two-dimensional models, often 
describing excitable media. 

For oscillatory systems, the motion in the concentration space largely takes place on a two- 
dimensional manifold. In particular, this condition applies near a supercritical Hopf bifurcation 
where, furthermore, the motion in the plane of oscillations can be described analytically to a 
very good approximation. Integrating the reaction-diffusion system close to a supercritical Hopf 
bifurcation can be a very time consuming and difficult task even on a large computer. For this 
reason, it is of great importance to introduce techniques that utilize the characteristics of the 
Hopf bifurcation in order to simplify the reaction-diffusion equation (RDE) without changing the 
essential properties of the solution. As shown by Kuramoto | 1984j| , such a description is offered 



by the complex Ginzburg-Landau equation (CGLE) which provides a systematic description of 
the chemical waves near a Hopf bifurcation (of the corresponding homogeneous system). He 
applied it to the Brusselator, a two-dimensional model for an abstract chemical system. 

The CGLE is an amplitude equation which describes the waves of the reaction-diffusion system 
in terms of a complex amplitude. It offers a number of advantages in addition to a reduction 
of the effective dimension of the original reaction-diffusion equation. It provides an intuitively 
natural description of oscillations which facilitates interpretations of results and, furthermore, 
introduces a scaling which makes it universal, depending only on two complex parameters but 
not e.g. on the distance from the bifurcation point. The scaling also solves a serious problem 
associated with the integration of reaction-diffusion equations, due to the critical slowing down of 
the motion in the concentration space which occurs near the bifurcation point. The description 
becomes increasingly precise as the Hopf bifurcation point is approached and is "exact" in the 
limit. 

As one moves further away from the bifurcation point the CGLE description gradually becomes 
less precise. Hence we wish to discuss how well the CGLE description approximates realistic 
reaction-diffusion systems representing actual chemical systems at finite distances from a Hopf 
bifurcation. This point is very important from an experimental point of view, since experimental 
observations of waves and spatio-temporal chaos require sufficiently large amplitudes and rates 
of development for spatial patterns. 
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The behavior of amplitude equations that apply close to a Hopf bifurcation have for example 
been studied by Huber et al. 1992 ] , Aranson et al. 1993 1, and Nicofis 1995| . However, to our 
knowledge, there have been no attempts to test how well the CGLE actually does approximate 
the behavior of a realistic reaction-diffusion system close to a supercritical Hopf bifurcation. By 
comparing results from simulations of a specific reaction-diffusion system with those of the asso- 
ciated CGLE, we show in this paper that the CGLE indeed gives a very satisfactory description 
of the original system. We focus attention on two specific conditions for an Oregonator model, 
where spiral patterns and chaotic patterns can occur when the corresponding homogeneous sys- 
tems exhibits simple oscillations with small amplitude. 



2 Review of Theory 

We consider a chemical reaction with n chemical components taking place under heterogeneous 
conditions described by the reaction-diffusion equation. The reaction-diffusion equation is de- 
fined by the rate expression and the diffusion matrix associated with the chemical system under 
consideration. Let the vector c describe the concentrations of the n dynamical species that par- 
ticipate in the reaction-diffusion system. The associated reaction-diffusion equation can then be 
written as 

c = f(c ;/ u)+D-V2 c , (1) 

where f is the rate expression and D is the diffusion matrix. Here /i designates some controllable 
parameter which can be varied experimentally and serves as a bifurcation parameter. Assume 
that this system has a homogeneous and stationary solution c s . The dynamics of small pertur- 
bations u(t, r; fj,) = c(t, r; fi) — c s (t; /x) of the stationary point can then be described through a 
Taylor expansion of the right-hand side of Eq. (|]) by 

u = J u + ^M:uu + ^Niuuu + . . . + D-V^u. (2) 

The matrix Jij = ^ is the Jacobian matrix, whereas M:uu and N:uuu denote quadratic and 
cubic forms in u, respectively defined by 



t 

j,k=l 



UjU k , (3a) 



iNluUU], = £ <?V ' 



UjiikUi. (3b) 



Cs 



M is often referred to as the Hessian. The stability of the stationary state c s to space-independent 
perturbations is determined by the eigenvalues of the Jacobian matrix associated with the lin- 
earization of the rate expression f(c;/x) in Eq. ([l|) around c s . In general, a bifurcation or loss of 
stability will occur if the real part of one or more of these eigenvalues change sign from negative 
to positive. Let us assume that c s looses stability via a supercritical Hopf bifurcation at \i = /xo 
due to a change of sign in the real part of a complex eigenvalue A = a + iu and its complex 
conjugate. At the bifurcation point, the eigenvalue A will therefore be purely imaginary, iloq. In 
the homogeneous system, this bifurcation gives rise to small sinusoidal oscillations of the con- 
centrations. In the limit fj, —*■ fJ,Q, the frequency of the oscillations will approach the value loq 
corresponding to the imaginary part of A at the bifurcation point. 
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The equations related to the Hopf bifurcation can be expanded in a formal parameter e which 
often is related to the amplitude of the oscillations. The quantities u, J(/i), M(/i), N(/z), 
A(/i) = a(fi) + iuj(fi) and the original parameter \x can be expanded as follows 

u = eui + e 2 u 2 + e 3 u 3 . . . , J = J + e 2 J 2 + e 4 J 4 . . . , 
M = M + e 2 M 2 + e 4 M 4 . . . , N = N + e 2 N 2 + e 4 N 4 . . . , 

A = A + e 2 A 2 + e 4 A 4 . . . , 1 ' 

\i = no + e 2 /i 2 + eV 4 . . . 

where Aj = aj + iujj,j = 0, 2, 4, . . .. The parameter e can be defined in many ways. In experimen- 
tal situations where the bifurcation parameter is a real physical quantity, it is often convenient to 
choose e in a dimensionless form and eliminate it from all expressions that are used to interpret 
results [Kosek et al. 1994 1. 



Consider the behavior of the system in states that locally are close to the homogeneous oscillatory 
solution. One can prove flKuramoto 1984 ] that close to the Hopf bifurcation point the dynamics 



of the n chemical species c(t, r; fj,) can be approximated by the following equation 

c(t, r;n)=c s + e (w(t, s)e iuJot U + + W(r, s)e~ iu ' ot TJ + ) . (5) 

Here U + is the complex right eigenvector associated with the bifurcating eigenvalue iwo at the 
bifurcation point e = (jx = fio)- The complex amplitude W(r, s) of Eq. (||) must satisfy the 
complex Ginzburg-Landau equation 

8W n n 

_ = \ 2 W - g\W\ 2 W + dV 2 s W, (6) 

which depends on t, r, and e through the time and space variables r and s, scaled with e as 
r = e 2 t and s = er. One purpose of the scaling is to obtain a description that is independent 
of the distance from the bifurcation point (as long as the approximation is applicable). The 
use of the scaled time r solves the problem of critical slowing down exhibited when the system 
approaches the bifurcation point. In other words, the closer we get to the bifurcation point, the 
longer the transient time for pattern development. As we shall see, this fact can indeed cause 
difficulties in numerical studies of the original reaction-diffusion equation, which are absent when 
the CGLE is used. 

The complex parameter d in the CGLE is determined by the equation 

d = U;-D-U+, (7) 

where U^j_ is the left eigenvector corresponding to the bifurcating eigenvalue i loq at the bifurcation 
point. The left eigenvector U?j_ has been normalized to fulfil the relation U^_-U + = 1. To find 
the complex coefficient g one first determines two vectors F 2 q and Fn by solving the two linear 



equations [[Kuramoto 1984|| 



(J -2iw I)-F 20 = -iM :U + U+, (8a) 
Jo-F n = -M :U + U+. (8b) 



The complex parameter g can be found by the expression | Kuramoto 1984 ] 



-U;-M :U + F n - U;-M :U + F 20 - ^U^ •N :U + U + U + . (9) 
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In Eq. (H) and Eq. (g), the tensors Jo, Mo, and No denote the values of J, M, and N at the 
bifurcation point fi = fiQ. 

The complex Ginzburg-Landau equation is often presented in the form 

W- (1 + if3)\W\ 2 W + (1 + ia)V 2 s W. (10) 



Or 

Here a = and = where d' = Red, = Imd, </ = Re g, and 5" = Img. One obtains 
this dimensionless version of the complex Ginzburg-Landau equation by performing the variable 
change r -> ct^V, s -> ^Js, W -> y|^exp(i ^r). 

One can easily show that in a one-dimensional spatial domain of infinite length Eq. ( |To| ) admits 
of plane wave solutions of the form 

W Q (x,t) = R Q e i ^ x ~^ t \ (11) 
where the frequency ujq and amplitude Rq have to satisfy the relations 



Rq = \Jl - Q 2 , 101 <1, (12a) 

= /?g 2 + a(i-g 2 ). (12b) 

By considering perturbations of the plane wave solutions, one can prove that these become 
unstable when a and (3 satisfy the relation 1 + a(3 < 0. The quantity 1 + a(3 = determines 
the boundary of the Benjamin-Feir instability which is associated with spatio-temporally chaotic 
behavior in the system. Since spiral wave solutions far from the spiral core asymptotically 
approach a plane wave locally, one would expect that crossing the boundary of the Benjamin- 
Feir instability also implies a breakdown of spiral waves and perhaps development of chaotic 
behavior. In Sec. |] we shall present examples showing the breakdown of spiral wave solutions in 
both the reaction-diffusion equation and the Ginzburg-Landau equation. 



3 Model description 

As a model for a real chemical system we have chosen the three dimensional Oregonator [ Field| 



fe Noyes 1974 ] . The Oregonator is based on the so-called FKN-mechanism, which provided the 



first successful explanation of the chemical oscillations that occur in the BZ-reaction [Field et al. 



1972 1 . During the last two decades, the Oregonator model has been modified in many ways by 
inclusion of additional chemical reaction steps or by changing the rate constants. The model 
used in this study consists of the original Oregonator mechanism 



(13a) 
(13b) 
(13c) 
(13d) 
(13e) 



BrO^ + Br" + 2H+ 




HBr0 2 + HBrO 


HBr0 2 + Br" 


fc 2; 


2HBrO 


2HBr0 2 


fc 3) 


HBr0 2 + BrOg + H 


BrO^ + HBr0 2 + H + 


&4 


2HBr0 2 + 2Ce 4+ 


Ce 4+ 


H } 


/Br" + Ce 3+ 



with rate constants k\, . 

factor. For notational simplicity we introduce A 



&5 from Field & Forsterling 1 1986 ] . The quantity / is a stoichiometric 

[Br0 3 - 



H 



[H+],X = [HBr0 2 ],y = [Br-], 
and Z = [Ce 4+ ], where [S] is the concentration of species S. For relevant experimental conditions 
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we may assume that [BrO^ ] and [H + ] are constant. By applying mass law kinetics, the time and 
space variation of the spatially distributed Oregonator model can be described by the following 
three coupled partial differential equations 



OX 

at 



— = k x AH 2 Y - k 2 HXY - 2k 3 X 2 + k A AHX + D x V 2 r X, 



8Y 
dt 



— = —k\AH 2 Y - k 2 HXY + k 5 fZ + D Y V 2 r Y, 



dZ 
dt 



— = 2k 4 AHX -k 5 Z + D z V 2 r Z, 



(14a) 



(14b) 



(14c) 



where Dx, Dy, and Dz are the diffusion constants of the species HBr02, Br~, and Ce + 
respectively (for dilute solutions, the diffusion matrix is diagonal to a good approximation). For 
a thorough discussion of the chemistry on which the Oregonator is based, the reader is referred 



to Tyson [1985]. The values of the stoichiometric factor / [Nielsen et al. 1991], of the rate 
constants and of the diffusion constants used in the numerical studies are shown in table [I]. For 



the diffusion constants, we have used estimates by Hynne & S0rensen [1993] 



4 Numerical simulations 

The parameters a and (5 in Eq. ( |i~0| ) have a complicated functional dependence on the parameters 
that characterize the Hopf bifurcation in the original reaction-diffusion system as is evident 
from Eqs. (|7|)-(||). If a parameter in the original model system is changed slightly, a new Hopf 
bifurcation can be found by adjusting another parameter, corresponding to the fact that the Hopf 
bifurcation is a generic co-dimension one bifurcation. The Ginzburg-Landau parameters a and 
P will also change due to Eqs. (0) and ©. It is therefore convenient to study the variation of the 
parameters a and along a branch in a Hopf bifurcation diagram where the Hopf bifurcation 
is continued as a function of two appropriate parameters in the particular system. For the 
present Oregonator model, we have chosen the concentrations of BrOg and H + as bifurcation 
parameters, since these easily can be varied experimentally and are almost constant during an 
experiment. The bifurcation diagram obtained for the Oregonator model is shown in Fig. |]. The 
curve consists of two branches representing super- and subcritical Hopf bifurcations shown as the 
solid and dashed parts of the curve in Fig. |]. This change of stability occurs at [BrOg ] = 0.269 M 
and [H + ] = 0.176 M (M = mol dm -3 ). 

Since Eq. (||) applies to a supercritical Hopf bifurcation only, our attention will be limited to 
the supercritical branch in the bifurcation diagram. One can calculate the parameters a and j3 
from Eqs. ([?]) and @ along the bifurcation curve plotted as functions of [Br0 3 ] in Fig. [| We 
see that the diffusion parameter a changes slightly as the Hopf bifurcation point varies with the 
bifurcation parameters [BrO^~] and [H + ]. This variation is due to the fact that the directions 
of the right and left eigenvectors of the Jacobian matrix Jo change with the Hopf bifurcation 
point. The coefficient (5 shows a more dramatic change along the bifurcation curve. At the 
right-hand edge of the supercritical branch, j3 diverges to infinity because the denominator of 
(3 = becomes zero at the transition from supercritical to subcritical behavior. An extended 
Ginzburg-Landau approach may still apply close to the singularity if a quintic term is added to 
the complex Ginzburg-Landau equation. However, this modification will not be considered in 
this paper. 

A crucial term which can be calculated from the values of a and (3 along the Hopf bifurca- 
tion branch is the Benjamin-Feir parameter 1 + a/3, which determines whether plane waves 
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are unstable (1 + af3 < 0) to small perturbations in one space dimension. The variation of 
1 + a/3 along the supercritical part of the Hopf bifurcation curve is exhibited in Fig. ||| which 
shows that a Benjamin- Feir instability starts at [BrOj ]bf = 0.073 M and [H+] B f = 0.348 M. For 
[BrOg ] < [BrOg^BF and [H + ] > [H + ]bf, we may therefore expect to find stable spiral wave so- 
lutions to the reaction-diffusion system as well as to the complex Ginzburg-Landau equation. 
Spiral waves are expected to become unstable when [BrOg] > [BrO^BF and [H + ] < [H + ]bf- 
(Recall that [BrOj] and [H + ] vary together on the Hopf bifurcation curve.) 

We have selected two points on the (supercritical) Hopf bifurcation curve for comparisons of the 
(Oregonator based) reaction-diffusion system with the corresponding complex Ginzburg-Landau 
equation as indicated on Fig. [|. The characteristic Ginzburg-Landau parameters for the two 
points are shown in table |2[ The points differ in the sign of the Benjamin-Feir parameter, 
1 + a/3; they are marked on Fig. ||. The complex Ginzburg-Landau equation applies near a Hopf 
bifurcation and is scaled so that its form is independent of the distance from the bifurcation 
point. Thus, a and (3 are defined at the bifurcation. However, we want to compare solutions to 
the complex Ginzburg-Landau equation which is independent of the distance from the bifurcation 
point with those of the reaction-diffusion system for which the distance from the bifurcation point 
is important. In the latter case, one must work at definite (finite) distances from the bifurcation 
curve. The points are chosen with the following considerations in mind: If a point is too close 
to the bifurcation curve, significant spatial changes in amplitude take so long (compared to a 
local period of oscillation) that the numeric solution of the reaction-diffusion equation becomes 
practically impossible. On the other hand, the points should be close enough to the bifurcation 
for the complex Ginzburg-Landau equation to apply to a reasonably good approximation. 

One way, by which one can estimate a reasonable choice for the distance from the bifurcation, 
is to study how well the Stuart-Landau equation, approximates the true limit cycle solution for 
the homogeneous chemical reaction system. The Stuart-Landau equation is simply the complex 
Ginzburg-Landau equation without the diffusion term, i.e. 

= \ 2 W-gW 2 W. (15) 

OT 

If we put W = Re 10 , then Eq. ( |l5|) has the following simple solution 



R = J-j, (16a) 




9(r) = {u 2 -o 2 9 —)t, (16b) 
9 

for the motion on the limit cycle. Expressing the solution in terms of the real time t we find that 
the Stuart-Landau prediction of the oscillations of the chemical concentrations will be described 
by the expression 

c(t) = c s + (e iL ^U+ + e-' lu ^ t V + ) , (17) 

where oj(fi) = ujq + (^2 — 0"2^t)M- By plotting and comparing Eq. (|l7]) together with the actual 
solution of the reaction system, one normally will get an initial estimation of whether or not the 
Ginzburg-Landau approach can be justified. 

As the first point of investigation, we have chosen [BrOg ] = 0.01205 M with [H + ] = 1 M, which 
lies on the stable side of the Benjamin-Feir boundary (1 + a(3 = 0.76). Plots of the Stuart-Landau 
solution and the actual solution to the homogeneous Oregonator model are presented in Fig. [|a. 
The point is chosen at a distance from the bifurcation where the amplitude of the Stuart-Landau 
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solution is 10% of the value of [Ce 4+ ] at the stationary point c s . The results for the reaction- 
diffusion equation and the CGLE on a two-dimensional spatial domain are shown in Fig. ||. One 
clearly sees that the spatial wavelengths of the spirals are in excellent agreement, 4.3195 cm and 
4.3194 cm respectively, determined by Fourier transformation. The amplitudes (not shown) also 
agree. Note also that the actual times it takes for the two spirals to evolve are almost the same, 
which clearly justifies the time scaling r = fit predicted by the Ginzburg-Landau theory. 

In order to see what happens when the distance from the bifurcation point is increased, we 
have chosen two different values of [BrOj] which give rise to amplitudes of 50% and 60% of the 
stationary concentration of Ce 4+ . The comparisons for the homogeneous systems, i.e. between 
the actual limit cycle and the Stuart-Landau prediction are shown in Figs. |]b and||]c. The results 
of the integration of the reaction-diffusion system and the CGLE for these two working points 
are presented in Figs. |6|a and |^b. From Fig. |6ja it is evident at such large value of the bifurcation 
parameter [BrOj] = 1.2950 x 10 -2 M, that the spiral center in the reaction-diffusion system no 
longer is localized in the center of the grid. This "symmetry break" is associated with the no- 
flux boundary conditions and we shall ignore it in our comparisons. Even though the complex 
Ginzburg-Landau equation fails in detail, it still reproduces the behavior of the reaction-diffusion 
system quite well. The spatial wavelengths are almost equal, 0.772 cm for the reaction-diffusion 
equation and 0.720 cm for the CGLE, and the characteristic time scales of spiral development 
are essentially the same for the two systems. 

We now consider the case [Br0 3 ] = 1.3596 x 10~ 2 M with an amplitude 60% of the average value 
for the homogeneous Stuart-Landau solution. From Fig. ^3 we clearly see that at this rather 
large distance from the bifurcation point, the Ginzburg-Landau approach breaks down. The 
CGLE description predicts a stable spiral solution for the reaction-diffusion system, which itself 
exhibits a totally different pattern: A small spiral tip is formed near the boundary. After a few 
windings the spiral becomes unable to maintain its structure. It becomes unstable and breaks 
up into a large number of small spiral cores distributed throughout the spatial domain of the 
system. The visual appearance of the solution to the reaction-diffusion equation does indeed 
look "complex". However, the system is not chaotic. This has been investigated by calculating 
the largest Lyapunov exponent A max which is zero. The system is thus characterized by an 
oscillatory state which in this case is dominated by a large collection of relatively small spirals. 
Even though the resemblance between the reaction-diffusion system and the CGLE seems to be 
lost completely for this particular situation, one should note that the spatial wavelength of the 
small reaction-diffusion spirals still are in good agreement with the spatial wavelength found in 
the single spiral wave solution to the CGLE. 

We now turn to a point on the bifurcation curve where the Benjamin- Feir parameter 1 + a/3 is 
negative, viz. [BrO^Hopf = 0.073 M and [H + ]nopf = 0.348 M where the Benjamin-Feir parameter 
equals 1 + a/3 = —0.13. The actual integrations of the reaction-diffusion system were carried out 
at [BrO^] = 0.0817 M with [H+] = [H + ] Ho pf- Here the Stuart-Landau amplitude is 10% of the 
average Ce 4+ concentration, see Fig. ||d. The results for the reaction-diffusion equation and the 
CGLE are shown in Figs. 0a (phase) and 0b (amplitude). In both simulations (RDE and CGLE) 
we see that initially a spiral wave begins to develop from the center of the region. However, when 
two spiral windings have been formed, the spiral becomes unstable and a shock-like circular wave 
develops from the spiral core. When this shock wave hits the boundaries of the spatial domain 
a number of new spiral cores are generated (these are often referred to as phaseless points). 
Some of these phaseless points will repeat the scenario just described, whereas others will be 
annihilated either at the boundaries of the system or by colliding with other phaseless points. 
Finally, we end up with a spatio-temporally chaotic state which is totally dominated by phaseless 
points. From Fig. 0a and we note the striking similarity between the behavior of the reaction- 
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diffusion equation and the corresponding CGLE. The development of the chaotic state takes 
place in almost the same time in the two systems. Furthermore, we see that the number of 
phaseless points in the spatial domain of the two systems are of the same order of magnitude. 
The chaotic character of the system is illustrated in Fig. ||, where the estimate of the largest 
Lyapunov exponent for the CGLE system, A max , is plotted as a function of time. For t — ► oo the 
curve converges to a limiting value, the actual value of the Lyapunov exponent. As is evident 
from Fig. ||, the Lyapunov exponent is positive corresponding to chaotic dynamics. 

The characteristic parameters used for the numerical integrations discussed in this section are 
summarized in table ||. 



5 Computational details 



The computation of the Hopf bifurcation diagram and the limit cycle solutions to the homo- 



geneous Oregonator model were done with the bifurcation analysis package CONT [Schreibei 
1995 1 . All integrations of the reaction-diffusion systems were performed on a CRAY92-computer 
using an explicit fourth order Runge-Kutta algorithm | Press et al. 1992 1 . The integration of 
the CGLE were done on a HPUX-workstation using the implicit Adams method from ODE- 



PACK | Hindmarsh 1982 1. All integrations were carried out with no-flux boundary conditions 
and initial conditions chosen as described in | Kuramoto 1984| , p. 106]. The Lyapunov exponent 



was computed by using a numerical approach described by Marek and Schreiber 1 109 1 ] 



6 Discussion 

We consider two ways of using the complex Ginzburg-Landau equation for chemical reaction 
diffusion systems. One is to simplify calculations with models. Another one is to base realistic 
modelling of actual reaction diffusion systems on the CGLE with parameters determined directly 
by experiments performed on homogeneous systems [Hynne & S0rensen, 1993]. 

The numerical simplification of the reaction-diffusion models by the CGLE is obtained in two 
ways. First, the CGLE effectively reduces the dimension of the state space to two, which saves 
computation time, particularly for high-dimensional models. For example, calculation of waves 
for the big model of Gyorgyi et al. [ 1990] with its 26-dimensional state space, is no more difficult 



numerically than for any other model — once the Ginzburg-Landau parameters have been cal- 
culated. (Here the main problem is to obtain all the relevant diffusion coefficients). Second and 
most importantly, the use of the CGLE solves the problem (with the RDE) of the critical slowing 
down of pattern development that occurs near a bifurcation. The CGLE works with amplitudes, 
and the major effect of the time oscillations is accounted for analytically (in the transforma- 
tion to the amplitude equation) rather than numerically. In fact, the critical slowing down has 
been completely scaled away in the CGLE. The importance of this feature can be appreciated 
by the significant reductions (by factors up to 1000) of the number of steps necessary for the 
integration, see table ||[ In practice this means that the CGLE-calculation can be made with 
a work station whereas a supercomputer may be needed for integrating the RDE. In addition, 
the CGLE provides the solution in the form of a time and space dependent complex amplitude. 
This form is very convenient for understanding solutions, and it can easily be transformed back 
to actual time and space dependent concentrations of the (multitude of) species participating in 
the reactions. 

The other way of using the CGLE for chemical systems is to set up the Ginzburg-Landau equation 
directly from experimental measurements on the corresponding homogeneous reaction system 
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using independently measured diffusion coefficients. It has been shown previously | Hynne & 



S0rensen 1993 | how the parameters can be obtained from quenching experiments |S0rensen & 



Hynne 198E ; Hynne et al. 1990| ; Vukojevic et al. 1993 ; Nagy et al. 1995 1. Thus, it is possible to 



model waves and chaos in real spatial systems without knowing the detailed mechanism of the 
chemical reactions. 

The previous discussion shows that the Ginzburg-Landau approach to chemical reaction-diffusion 
systems is potentially very important. However, the use of the CGLE depends on the validity of 
the approximation for any particular system and operating point. 

To throw some light on the question of the validity of the complex Ginzburg-Landau approx- 
imation we have compared solutions to a model reaction-diffusion system with their CGLE- 
approximations. For a chemical reaction diffusion system based on the Oregonator model of 
the Belousov-Zhabotinsky reaction, we find that there exist substantial regions for which the 
CGLE provides an excellent description of the dynamical behavior. This conclusion applies to 
parameter regions where the solutions show spatio-temporal chaos as well as regions with stable 
spiral waves. In the latter case, we demonstrate how the CGLE approximation eventually fails 
beyond a certain limit. 

Experimentally, it is difficult to work too close to a Hopf bifurcation because small amplitude 
waves are difficult to detect and because of slow pattern evolution (critical slowing down). We 
conclude from the present study that one can reasonably expect a Ginzburg-Landau descrip- 
tion to remain applicable to a system at experimentally feasible distances from the bifurcation 
point, provided the properties of the Oregonator model can be regarded as representative of the 
investigated system. 

One condition for using the conclusions of the present study is that the real parts of the eigen- 
values of the Jacobian matrix (other than the bifurcating pair) are sufficiently negative as in 
the Oregonator model considered here so that excitations out of the plane of oscillations decay 
fast. We are currently studying what happens when one such "transient mode" is slow, and the 
present work serves as a preliminary to that study as well as a prerequisite to using the CGLE 
with experiments and models. 

To conclude, waves and spatio-temporal chaos in real chemical systems can be described by a 
complex Ginzburg-Landau equation with parameters a and j3 determined by quenching experi- 
ments performed on the specific reaction at the specific point of operation used. The paper has 
demonstrated that, typically, the Ginzburg-Landau description will apply in a sufficiently wide 
region surrounding the bifurcation set to make the description useful in practice. 
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Captions of tables and figures 

Table [l]: Rate constants ki,...,ks, stoichiometric factor / and diffusion constants Dx, Dy and 
Dz used in the numerical integration of the Oregonator model, Eq. (|l4|). The values of k%, . . . , fej 
are from Field and Forsterling [[1986(1 , k§ and / from Nielsen et al. jl991| and diffusion constants 
from Hynne Sz S0rensen | 1993| , (M = mol dm -3 ). 



Table |2|: Parameter values of the complex Ginzburg-Landau equation at the two points of 
the (1 + a/3)-curve plotted in Fig. ||| corresponding to the values of the two-parameter point 
([BrOg^Hopf j [H + ]Hopf)- The values of [BrO^] actually used in the numerical integration of the 
reaction-diffusion system are given in table ||[ Note that these values are off the bifurcation curve. 
The number ujq is the frequency of the sinusoidal oscillations exactly at the Hopf bifurcation 
point, whereas &i and 0J2 are the derivatives of the real and imaginary parts of the eigenvalue 
A with respect to the bifurcation parameter (i.e. [BrO^] for this particular example). The 
parameters g and d are the coefficients of the nonlinear term and the diffusion term in the complex 
Ginzburg-Landau equation, from which the dimensionless parameters a and (3 are calculated (see 
Eqs. (^) and (10)). The Benjamin-Feir parameter 1 + af3 for the two Hopf bifurcation points 



have been selected in order to study two different situations: At the first point (1 + aji > 0) we 
find that the reaction-diffusion system to admits of stable spiral wave solutions, whereas for the 
second point (1 + af3 < 0) we find chaotic behavior. 

Table Integration parameters for the reaction-diffusion equation and the corresponding com- 
plex Ginzburg-Landau equation for different amplitudes relative to the stationary Ce 4+ concen- 
tration and for the two selected points on the Hopf bifurcation curve (Fig. [I]). The table shows 
the final integration time, t e nd ; an d size, (r x ,r y ), of the 2-dimensional domain for the specific bi- 
furcation point distances of [BrOg ] used in the integration of the reaction-diffusion system. r en d 
and (s x ,s y ) indicate the corresponding scaled values used in the integration of the CGLE. The 
size of the grid is also shown together with a comparison between the steps used by the solver 
to reach t en d for the reaction-diffusion system and the associated CGLE (equal grid sizes was 
chosen for the integration of a specific reaction-diffusion system and the corresponding CGLE). 

Figure [l]: Bifurcation diagram showing the locations of Hopf bifurcations in the Oregonator 
model in the plane of the two parameters [BrOg ] and [H + ]. At a specific point [BrOg ] = 0.269 M 
and [H + ] = 0.176 M, the Hopf bifurcation changes stability from supercritical to subcritical corre- 
sponding to the solid and dashed curves respectively. For each point on the solid curve, a unique 
Ginzburg-Landau equation is defined, since the Ginzburg-Landau parameters ui2, 02, 9 and d 
can be calculated from Eqs. (0)-®. The comparisons between the reaction-diffusion system and 
the complex Ginzburg-Landau equation have been carried out at the two points indicated on the 
curve. 

Figure |||: (a): The dimensionless diffusion parameter a along the supercritical branch in Fig. ||, 
plotted as a function of the parameter [BrOcT]. (b): Variation of the dimensionless parame- 
ter P associated with the nonlinear term in the complex Ginzburg-Landau equation along the 
supercritical branch in Fig. [j] plotted function of the parameter [Br0 3 ]. 

Figure ^: The figure shows the variation of the Benjamin-Feir stability parameter 1 + a(3 along 
the curve of Hopf bifurcation points in the Oregonator model with [BrOj ] (a) and [H + ] (b) 
as parameter. From (a) and (b) we see that the Oregonator reaction-diffusion system crosses 
the Benjamin-Feir instability border at the point [BrO^~] = 0.073 M and [H + ] = 0.348 M where 
1 + af3 changes sign. As reference for further numerical investigations we have chosen two points 
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indicated on figure (a) and (b) with a filled and an open circle: One point showing stable spiral 
wave solutions and one point exhibiting spatio-temporally chaotic behavior. 

Figure Comparisons of the actual periodic solutions to the homogeneous Oregonator model 
(solid lines) and its approximation by the Stuart-Landau equation (dashed lines) at various 
distances from the bifurcation point, (a)-(c): Periodic solutions with amplitudes of 10%, 50% and 
60% of the stationary concentration of Ce 4+ at the Hopf bifurcation point represented in Fig. [l] 
with a filled circle, (d): Oscillations with an amplitude of 10% of the stationary concentration 
of Ce 4+ associated with the Hopf bifurcation point represented with an open circle in Fig. |]. 

Figure |5|: The phase (ft of the oscillations for a stable spiral solution to the Oregonator reaction- 
diffusion system (RDE) and the associated complex Ginzburg-Landau equation (CGLE). The 
parameters used in the integrations correspond to the point marked with a filled circle in 
Figs, [i] and [3|. In this case there is a very good agreement between the time scale under which 
the two patterns develop as well as the spatial wavelength of the two spirals. The numbers 
below each of the six snapshots of the RDE shows the time elapsed for that particular state 
in the RDE in units of 10 5 s. The corresponding snapshots for the CGLE are made at corre- 
sponding scaled times. The reaction-diffusion system has been solved for the parameter value 
[BrOg ] = 1.2051 x 10~ 2 M corresponding to an amplitude of 10% of the stationary concentration 
of Ce 4+ for the homogeneous Oregonator system. 

Figure Hf Deviations between the reaction-diffusion system and the CGLE description become 
more evident when a larger distance from the bifurcation point is chosen. The figure compares 
the two systems at two finite distances corresponding to an amplitude of 50% (a) and 60% (b) 
of the stationary Ce 4+ concentration. Due to nonlinear effects which become dominating as the 
distance increases the spiral moves away from the center of the domain and freezes at a point 
close to the boundary of the domain. In (a) the general characteristics (spatial wavelength and 
time of evolution) of the spirals in the two systems are still in good agreement. When the distance 
is increased further (b) the resemblance between the reaction-diffusion system and the CGLE is 
completely lost. The final state of the reaction-diffusion system is characterized by a number of 
small spirals distributed throughout the domain of the system. However, the spatial wavelength 
of these spirals are still described well by the CGLE. 

Figure |7|: Development of a chaotic state in the Oregonator reaction-diffusion system (RDE) 
and the associated complex Ginzburg-Landau equation (CGLE) represented by the phase (ft (a) 
and the amplitude (b) of the oscillations of [Ce 4+ ]. The characteristic parameter values used 
in the RDE and the CGLE correspond to the parameter point marked with an open circle in 
Figs. [I] and || where the Benjamin-Feir parameter satisfies 1 + a[3 < 0. Indications below the 
twelve snapshots of the RDE shows the evolution time for that particular state in the RDE as 
well as for the equivalent state in the CGLE measured in units of 10 5 s. Initially a spiral is 
formed at the center of the grid, but after approximately two windings have been, formed this 
structure breaks down being replaced by a set of phaseless points. (A phaseless point appears 
as a very localized dark spot on the [Ce 4+ ]-plot). The phaseless points will either split into 
other phaseless points or be annihilated at the boundary of the spatial domain. This scenario is 
then repeated, resulting in a very complicated interaction between the phaseless points. Finally 
a spatio-temporally chaotic state is reached in both the RDE and the CGLE. In the plots 
of the phase and the amplitude of [Ce 4+ ] we see a striking similarity between the behavior 
exhibited by the RDE and the CGLE. The two patterns develop in an almost synchronized 
manner. Even when the spiral is fully destroyed (t = 2.5 x 10 5 s) the positions of the phaseless 
points are still very similar in the corresponding RDE and CGLE states. The synchronization 
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is however lost at the final chaotic state, due to the fact that this chaotic state is characterized 
by a critical dependence on the initial conditions. The reaction-diffusion system has been solved 
for the parameter value [BrOcT] = 8.1935 x 1CT 2 M corresponding to an amplitude of 10% of the 
stationary Ce 4 + concentration. 

Figure H Largest Lyapunov exponent, A max , calculated for the complex Ginzburg-Landau 
equation for the state presented in Fig. ^ (shown in Figs. |l| and || with an open circle). The 
value of A max is plotted as a function of time for a point where 1 + af3 < giving rise to spatio- 
temporal chaos and a positive Lyapunov exponent. In this sense, the state of the system can be 
characterized as chaotic. 
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